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ABSTRACT 

Popular cosmological scenarios predict that galaxies form hierarchically from the merger of many progeni- 
tors, each with their own unique star formation history (SFH). We use a sophisticated approach to constrain the 
SFHs of 4517 blue (presumably star-forming) galaxies with spectroscopic redshifts in the range 0.2 < z < 1.4 
from the All- Wavelength Extended Groth Strip International Survey (AEGIS). This consists in the Bayesian 
analysis of the observed galaxy spectral energy distributions with a comprehensive library of synthetic spectra 
assembled using realistic, hierarchical star formation and chemical enrichment histories from cosmological 
simulations. We constrain the SFH of each galaxy in our sample by comparing the observed fluxes in the B, 
R, I and Kg bands and rest-frame optical emission-line luminosities with those of one million model spectral 
energy distributions. We explore the dependence of the resulting SFHs on galaxy stellar mass and redshift. 
We find that the average SFHs of high-mass galaxies rise and fall in a roughly symmetric bell-shaped manner, 
while those of low-mass galaxies rise progressively in time, consistent with the typically stronger activity of 
star formation in low-mass compared to high-mass galaxies. For galaxies of all masses, the star formation 
activity rises more rapidly at high than at low redshift. These findings imply that the standard approximation of 
exponentially declining SFHs widely used to interpret observed galaxy spectral energy distributions may not 
be appropriate to constrain the physical parameters of star-forming galaxies at intermediate redshifts. 
Subject headings: galaxies: evolution - galaxies: star formation - galaxies: stellar content 



1. INTRODUCTION 

Constraints on the stellar content of galaxies are often de- 
rived from multi wavelength observations by assuming that the 
star formation history (SFH) of an individual galaxy can be 
approximated by a simple declining exponential function of 
the form exp(— f/r), where t is the galaxy age and r the star 
formation timescale. Such "r-models" have been used suc- 
cessfully to estimate, for example, the stellar masses of nearby 
spiral galax ies from fits of rest- frame optical/near-infrared 
colors (e.g., iBell & de Jong 20Mh . Despite this success, an 
increasing number of analyses have pointed out the limita- 
tions of T-models, particularly for applications to studies of 
high-redshift galaxies. T his is the case of Lym an-break galax- 
ies at redshifts z > 2 jPapovich et al. 20011) an d blue star- 
forming galaxies at z ~ 2 (IShapley et al. 2005h . for which 
the addition of another stellar component on top of a r-model 
can change the derived stellar masses by a factor of several. In 
fact, the major episodes of star formation in ultraviolet (UV)- 
luminous Ly man-break galaxies could last only a few hundred 
milUon years dStark et al. 20091) . Reddy et al. (2012) also find 
that, for galaxies at 1.5 < z < 4, fits of broadband spec- 
tral energy distributions (SEDs) using exponentially declin- 
ing SFHs produce on average star formation rates (SFRs) 5 to 
10 times lower than those inferred from more direct estimates 
based on the combined analysis of U V and mid-infrared emis- 
sion. They conclude that SFHs in which the SFR rises with 
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time are more appropriate than declining r-models for these 
galaxies. This is reinforced by the conclusions of Maras- 
ton et al. (2010) and Pforr et al. (2012) that SED fits of 
galaxies at z ^ 2 using declining r-models produce unre- 
alistically low ages, and that exponentially rising r-models 
should be preferred for these galaxies. In fact, simulations 
of galaxy formation in a hierarchical universe predict com- 
plex, stochastic SFHs of individual galaxies, which are not 
well approximated by simple decli ning exponential laws (e.g., 
iLee et al. 20091: IWuvts et al. 20091) . 

To better characterize the SFHs of galaxies at moderate 
and high redshift, we require more sophisticated, physically 
motivated spectral analysis tools. In this Letter, we achieve 
this goal by appealing to the approach recently proposed by 
Pacifici et al. (2012) to constrain the SFHs of 4517 galaxies 
with secure spectroscopic redshifts in the range 0.2 < z < 
1.4. Specifically, we fit the observed broadband SEDs and 
emission-line strengths of these galaxies using a comprehen- 
sive library of synthetic spectra. We assemble this library by 
combining state-of-the-art models of stellar population syn- 
thesis, nebular emission, and attenuation by dust with star 
formation and chemical enrichment histories derived from 
the semi-analytic post-treatment of a large-scale cosmological 
simulation. This approach allows us to characterize the SFHs 
of galaxies in different mass and redshift ranges in a physi- 
cally motivated and statistically reliable manner, using realis- 
tic, hierarchical mass-assembly histories. The resulting best- 
estimate SFHs differ significantly from simple exponentially 
declining laws, and the dispersion about the average SFH in 
any mass and redshift range can be quite large. 

The paper is organized as follows. In |J2 we present ob- 
servations of the B-, R-, I- and A'^-band fluxes and rest-frame 
optical emission-line luminosities of 4517 blue galaxies from 
the All- Wave length Extended G roth Strip International Sur- 
vey (AEGIS, iDavis et al. 2007h . In we describe the li- 
brary of synthetic SEDs built to interpret these observations. 



2 



Pacifici et al. 



We report our results in S|4]and present our conclusions in [jS] 
Throughout this Letter, we adopt a standard ACDM cosmol- 
ogy with JIm = 0.3 , r?A = 0-7, h = 0.7, and a Chabrier 
initial mass function (IChabrier 2003h . 

2. DATA 

To explore the SFHs of galaxies in wide ranges of stel- 
lar mass and SFR at different cosmic epochs, we appeal to 
the AEGIS survey. This combines photometric and spec- 
troscopic observations of the Extended Groth Strip at wave- 
lengths over the whole range from X-ray to radio. In this 
initial analysis, we restrict ourselves to optical and infrared 
wavelengths. This is because multiband observations at wave- 
lengths spanning from blue-ward the 4000A break out to the 
rest-frame /-band set valuable constraints on a galaxy's SFH, 
while spectroscopically measured emission lines are required 
for s olid redshift determinati ons and precise SFR estimates 
(e.g.. lBrinchmann et al. 20041) . From the AEGIS catalog, we 
extract 6246 gal axies wit h available photom etry at B, R, I 
dCoil et al. 20041) and dBundv et al. 2006h . The resulting 
sample is magnitude-limited to ^"^=22. 5. From this sample, 
we select 4517 potentially star-forming galaxies {U — B < 
1.0), which do not show contamination by an Active Galactic 
Nucleus (AGN). In particular, we reject the sources detected 
by Chandra in addition to those identified as AGNs follow- 
ing the criteria by Kauffmann et al. (2003) (rest-frame opti- 
cal emission-line ratios), Yan et al. (20 11) (rest-frame optical 
emission-line ratio versus color), and Donley et al. (2007) 
(slope of the Spitzer IRAC SED). Since only 2% of the poten- 
tially star-forming galaxies are detected by Chandra, remov- 
ing this small percentage does not bias our results. For all 
4517 galaxies, Keck/DEIMOS spectra a re avai lable from the 
DEEP2 Redshift Survey (Newman et al . 20121) . They cover 
the wavelength range 6500-9 lOOA at a resolution of 1.4A full 
width at half-maximum (FWHM). The accuracy of the red- 
shifts is 30 km s^^. Such resolution allows us also to extract 
reliable fluxes of [Oii1AA3 726.G, 3728.8, H/3, [Oiii]A5007, 
and Ha dWeiner et al. 20071) . when these lines fall into the ob- 
served wavelength range of the spectra (Table [T]!. This set of 
photometric and spectroscopic data probes a rest-frame wave- 
length range from A ~ 3300A-1.8/xm at z = 0.2 to A 
1650-9150Aatz = 1.4. 

3. MODELING APPROACH 

We use the recently developed spectral analysis tool of 
Pacifici et al. (2012) to derive optimal constraints on the SFHs 
and stellar masses of the galaxies in the sample described 
above. This tool relies on the combination of a large library 
of physically motivated SFHs from cosmological simulations, 
with state-of-art models of stellar population synthesis, nebu- 
lar emission, and attenuation by dust. 

3.1. Library of model spectral energy distributions 

Following Pacifici et al. (2012), we build a comprehen- 
sive library of model star formation and chemical enrich- 
ment histories by performing a post-treatment of th e Mil- 
lennium cosmological simulation ( iSpringel et al. 20 05^ using 
the semi-analytic models of De Lucia & Blaizot (2007). For 
the present analysis of 4 broadband photometric fluxes and 
a few emission lines, we draw only 100,000 galaxy SFHs 
from z = 127 to the present time. By analogy with Paci- 
fici et al. (2012), we broaden the parameter space by redraw- 
ing, for each selected SFH, 10 different combinations of the 



following parameters: redshift of observation (in the range 
0.1 < z < 1.5); evolutio nary stage (see Sections 2.1 and 
3.1.2 in lPacifici et al. 20"T2|) : current (i.e., averaged over a pe- 
riod of 10 Myr before the galaxy is observed) specific SFR 
(—2 < log(V's/Gj/r~^) < 1); and current gas-phase oxygen 
abundance (7 < 12 + log(0/H) < 9.4). 

We generate a library of 1 million galaxy SEDs by com- 
bining this library of star formation and chemical enrichment 
histories with the latest version of the Bruzual & Chariot 
(2003) stellar population synthesis models, the galaxy neb- 
ular emission model of Chariot & L onghetti (20 01) (based on 
the photoionization code CLOUDY: iFerland 199 6). and the 2- 
component dust model of Chariot & Fall (2000). We take the 
same probability distribution of the total optical depth of the 
dust (fy) as in Pacifici et al. (2012), but with an upper limit 
of fv =3 (instead of 4) which is more appropriate for the 
current sample of blue star-forming galaxies. We also draw 
randomly the slope of the attenuation curve in the interstellar 
medium in the range 0.4-1 . 1 . 

This procedure allows us to go beyond the typical approach 
of using simple idealized functions at fixed metallicity to de- 
scribe the SFHs of galaxies. A potential drawback is that 
the star formation and chemical enrichment prescriptions are 
drawn from a specific semi-analytic model. Our resampling 
of the library, as described above, minimizes the impact of 
this choice. 

3.2. Fitting procedure 

We use the Bayesian approach detailed in Pacifici et al. 
(2012) to constrain the stellar masses, SFRs and SFHs of 
the 4517 galaxies in the observational sample of Sj2]using the 
above library of one million model SEDs. For each galaxy in 
the sample, we compute the likelihood that each model in the 
library reproduces the observed SED. Specifically, we com- 
pare the observed and modelled photometric fluxes in the B, 
R, I and Kg bands (observer frame) and the emission-line lu- 
minosities of [O II], H/3, [O III] and/or Ha (for each flux mea- 
surement, we adopt the errors from the AEGIS catalog and 
impose a minimum error of 10%). The exact lines fitted de- 
pend on the rest-frame wavelength coverage of the spectra, 
as listed in Table [T] Since the observed galaxies have ac- 
curate spectroscopic redshifts, we compute likelihoods only 
for the model galaxies lying within 0.05 of the spectroscopic 
redshift. In practice, this amounts to fitting each observed 
galaxy with ^ 71, 500 models. For each observed galaxy, we 
use the computed model likelihoods to build the probability 
density functions of stellar mass, il/^,, and SFR, ip (normal- 
ized to the absolute observed /iTs-band luminosity). We take 
the best estimates of M,^ and i/j to be the median values of 
these distributions. The typical uncertainty (defined as half 
the 16th-84th percentile range of the cumulative probability 
density function) is around 0.1 for Mi, and 0.3 for ijj. Actual 
average uncertainties in bins of stellar mass and redshift are 
indicated in parentheses in Table [T] For each galaxy in the 
sample, we also retain the best-fit SFH. 

4. STAR FORMATION HISTORIES OF BLUE GALAXIES 

Figure [T] illustrates the constraints derived on the star for- 
mation histories of blue galaxies at redshift around unity 
using the approach described in f|3] We show the best 
estimates of -0 as a function of M,^ for the 900 galax- 
ies at redshifts between 0.85 and 1.1 in the sample of ij2] 
(grey points). The results are consistent with the stan- 
dard "star-formation main sequence," which describes the 
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TABLE 1 

Median stellar masses and SFRs of AEGIS galaxies in different redshift bins (average uncertainties on individual 

MEASUREMENTS IN EACH BIN ARE INDICATED IN PARENTHESES). 



Redshift range 


Observed emission lines 


log(M*/M0) 


log[V'/(AfQyr-i)l 


Number of galaxies 






Low-mass" 


High-mass 


Low-mass" 


High-mass 


Low-mass" 


High-mass 


0.20 <z< 0.45 


H/3, [O III], Ha 


9.37(0.11) 


10.48(0.13) 


-0.185(0.308) 


0.375(0.347) 


167 


16 


0.45 <z< 0.70 


H/3, [0 III] 


9.40(0.09) 


10.56(0.12) 


-0.065(0.266) 


0.765(0.387) 


272 


68 


0.70 <z< 0.85 


[Oii],H/3, [Oiii] 


9.41(0.08) 


10.54(0.11) 


0.115(0.240) 


0.855(0.396) 


261 


124 


0.85 <z< 1.10 


[Oii],H/3 


9.43(0.09) 


10.52(0.11) 


0.625(0.301) 


1.155(0.546) 


191 


105 


1.10 <z< 1.40 


[On] 


9.50(0.10) 


10.54(0.10) 


0.965(0.229) 


1.335(0.477) 


96 


121 



" Galaxies with stellar masses between 1.6 and 4.0 x 10" Af© (blue points in Figure[2](. 
Galaxies with stellar masses between 2.5 and 6.3 x 10^" A/q (green points in Figure|2](. 
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Fig. 1. — Best estimates of SFR and stellar mass for the 900 galaxies at 
redshifts between 0.85 and 1.1 in the sample of [|2](grey points). The blue 
dash-dotted line represents the best-fit model SFH of an example galaxy of 
low stellar mass, while the green solid line represents the best-fit model SFH 
of an example galaxy of high stellar mass. Inset panels show SFR versus 
lookback time for these two example galaxies (at each time step, the SFR in- 
cludes contributions from all components of the merger tree that end up in the 
galaxy by the redshift of observation). Both best-fit SFHs display significant 
fluctuations in time, and neither would be well approximated by a r-model. 



systematic incre ase of the SFR with s t ellar mass and red- 



shift (see e.g., |Brinchmann et al. 2004 


Noeske et al. 20071 


iDaddi et al. 20071 IWhitaker et al. 2012bl ' 


Also shown in the 



inset panels in Figure [T] are the best-fit SFHs of two ex- 
ample galaxies at similar redshift, z « 0.9, with similar 
current SFR, i}] k, 5M0yr~^, but different stellar masses, 
= 4.3 X IO^Mq and 2.9 x W°Mq. It is interesting 
to note that the evolutionary paths of these galaxies remain 
within the star-formation main sequence over their entire life- 
times. The two SFHs show the characteristic fluctuations aris- 
ing from the merger histories of the dark matter halos and the 

' The normalization of the star-formation main sequence in Figure [T] is 
about 0.2 dex lower than that found by Noeske et al. (2007). This can be 
attributed to the different prescriptions adopted to derive SFR estimates and 
account for dust attenuation in that study. 



semi-analytic prescriptions used to model the baryons (gas in- 
fall, cooling, star formation, feedback). Neither SFH would 
be well approximated by a r-model. 

We now explore the average SFHs of galaxies in different 
bins of stellar mass and redshift. Figure |2] show the star- 
formation main sequence in five redshift bins spanning the 
range 0.2 < z < 1.4. Along these relations, we select two 
stellar-mass bins sampled at all redshifts: a low-mass bin 
(1.6-4.0 X IO^Mq) and a high-mass bin (2.5-6.3 x 10^° Mq). 
Table[T]lists the median stellar masses and SFRs in these two 
bins at all redshifts, together with the number of galaxies in 
each bin. For each redshift range. Figures |3] and |4] show the 
average best-fit SFHs of the galaxies in the low- and high- 
mass bins, respectively, as well as the root mean square de- 
viations about these averages. In each case, the average SFH 
is computed by first recording the best-fit SFH of each galaxy 
at cosmic times prior to the redshift of observation, and then 
co-adding the individual SFHs of the different galaxies. The 
resulting SFH is normalized at each time step by dividing by 
the number of galaxies forming stars at that time. The wig- 
gles in the SFHs of individual galaxies (see Figure [T]i tend to 
be smoothed out by such averages, except in bins containing 
few galaxies (left-most panel in Figure |4]i. We find that vir- 
tually identical average SFHs are obtained when selecting the 
second or third best-fit SFH for each galaxy. 

The SFR of low-mass galaxies (FigureO increases steadily 
from early times to the epoch of observation, while that of 
high-mass galaxies (Figure|4]i first rises and then declines in a 
bell-shaped manner Remarkably, at a fixed stellar mass, the 
shape of the average SFH is independent of galaxy redshift, 
although the evolution proceeds faster in the most distant 
galaxies. The lookback time at which galaxies had formed 
half their current stellar mass at a fixed redshift (vertical 
dashed lines in Figures |3] and |4]i increases with galaxy mass. 
This implies that the stars in massive galaxies in our sample 
are always typically older than those in low-mass galaxies. 
We note that the model library described in SjTT] includes a 
comprehensive range of SFH shapes from declining, to rising, 
to bell-shaped, to erratically bursty at any redshift. Hence the 
results of Figures |3] and |4] do not simply reflect the prior dis- 
tribution of SFHs in the model library. In fact, as shown by 
Figures|3]and|4]themselves, there is a significant dispersion in 
the SFHs of individual galaxies. 

It is of interest to note that the average SFHs of the galax- 
ies in our sample can be approximated by "delayed r-models" 
dBruzual & Kron 19801) of the type exp (—t/r), with t span- 
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Fig. 2. — Star-formation main sequence (i.e., SFR versus stellar mass) for the galaxies in five redshift bins spanning the range 0.2 < 2 < 1.4 in the sample of 
^ A bin of low-mass galaxies (M* = 1.6-4.0 X 10^ Mq) is highlighted in blue, and one of high-mass galaxies (M* = 2.5-6.3 X 10^" Mq) highlighted in 
green. 




lookback time / Gyr 



Fig. 3. — Average best-fit SFHs of the galaxies in the low-mass bin of Figure|2] in the same redshift bins spanning the range 0.2 < z < 1.4 from left to right 
(black solid line), and root mean square deviations about these averages (blue shadow). At any redshift, the average SFH of low-mass galaxies keeps rising in an 
extended way, inconsistent with declining r-models. 




lookback time / Gyr 



Fig. 4. — Average best-fit SFHs of the galaxies in the high-mass bin of Figure[2] in the same redshift bins spanning the range 0.2 < z < 1.4 from left to 
right (black solid line), and root mean square deviations about these averages (green shadow). At any redshift, the average SFH of high-mass galaxies rises and 
falls in a bell-shaped manner, which cannot be approximated by simple exponentially declining nor rising r-models. These findings imply that, at any redshift, 
low-mass galaxies (Figure[3j have on average higher cun'ent specific SFR (i/jg = ip/Mi,) than high-mass galaxies. 
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ning the range from 0.2 to 1.0 (from high to low redshift) for 
the high-mass galaxies and from 0.4 to 3.0 for the low-mass 
ones. The exponent d defines the ratio of the time of maxi- 
mum SFR to the star formation timescale {d = tj^ax/T). Of 
course, fluctuations arising from the specific merger history 
of any individual galaxy will make it deviate from the average 
SFH (Figure[B. 

From the best-fit SFHs extracted for each galaxy in our 
analysis, we can also make predictions for the star-formation 
main sequence at earlier times. We find that at z w 2 the 
predicted dispersion about the median SFR in different bins 
of stellar mass along the sequence is roughly consistent with 
observations by e.g., Daddi et al. (2007), Reddy et al. (2012), 
and Whi taker et al. (2012). As expected from the star- 
formation main sequence we observe at z < 1.4 (see foot- 
note |7]i, the normalization of the sequence is sUghtly lower 
than in these other studies. 

5. SUMMARY AND CONCLUSION 

We have measured the SFHs of 4517 blue galaxies at red- 
shifts between 0.2 and 1.4 in the AEGIS Survey using a new, 
sophisticated spectral analysis tool. To achieve this, we built a 
comprehensive library of galaxy SEDs by combining a semi- 
analytic post^ti^atjnentofthe Millennium cosmological sim- 
ulation (ISpringel et al. 20051) with state-of-the-art models of 
stellar population synthesis, nebular emission, and attenua- 
tion by dust. For each galaxy in our sample, we retrieved the 
best-fit model SFH and the probability density functions of 
stellar mass and SFR using a Bayesian analysis. We reach the 
following conclusions: 

• The best-fit SFHs of individual galaxies exhibit char- 
acteristic fluctuations arising from the merger histories, 
which cannot be well described by simple r-models. 

• At all redshifts in the explored range 0.2 < z < 1.4, 
the average SFH of high-mass blue galaxies (^ 4 x 
IO^^A/q) rises and falls in a bell-shaped manner, sug- 
gesting that these galaxies are gradually turning off 
their star formation activity. This behavior cannot be 
approximated by simple exponentially declining nor 
rising r-models. 

• The average SFH of low-mass galaxies (- SxlO^M©) 
keeps rising in a more extended way, also inconsistent 
with declining r-models. 

• The average SFHs of both high- and low-mass galax- 



ies can be reasonably well approximated by delayed r- 
models of the form f'* exp(— i/r). 

• The shapes of the SFHs suggest that, at any redshift, 
low-mass galaxies have on average higher specific SFR 
(tjjg = i/j/Mi,) than high-mass galaxies, consistent 
with the trends in the s tar-formation main sequence 
(e.g., iNoeske et al. 20071) and in the latest simulations 
ofBehroozietal. (2012). 

• At fixed stellar mass, the shape of the average SFH is 
independent of galaxy redshift, although the evolution 
proceeds faster in the most distant galaxies. 

• At fixed redshift, the lookback time by which half the 
stellar mass of a galaxy has formed (vertical dashed 
lines in Figures [3] and |4]i increases with stellar mass. 

• The SFHs retrieved from our analysis tend to make 
the galaxies remain within the star-formation main se- 
quence for their entire lifetimes. 

The results presented in this Letter will be extended in the 
near future by exploiting optical stellar absorption-line signa- 
tures and infrared photometry to improve the constraints on 
the SFHs (Pacifici et al. in preparation). 
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